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SECTION 1 


INTRODUCTION 


1.1 Simulation Extensions 

The objective of this project is to identify and quantify those generic 
mission, scenario and sensor parameters (as well as their interactions) which 
drive the performance of space-based staring IR surveillance systems. These 
efforts have been divided into two efforts: Simulation Enhancements; and 
Simulation Applications. 

Work in the Simulation Enhancements area is directed toward general¬ 
izing those models in the Draper Integrated Simulations (DIS) which are too 
limited in scope or perhaps restricted to one particular system. Models have 
been prepared for the platform, sensor (including focal plane and signal pro¬ 
cessor), scene, and techniques for interfacing the DIS with an in-house image 
processor for graphics display. 

The Simulation Applications effort studies issues related to generic 
surveillance system performance drivers using tools from the Simulation 
Enhancement effort. Typical examples are 3ensor line-of-sight stability, 
focal D*, and critical scene/scenario characteristics. 


HALO Optics 


HALO is a multi-mirror large optical imaging system with a wide field 
of view. To maintain alignment as well as image quality, the mirrors are 
actively controlled with actuators. The term "deconvolution" implies a sens¬ 
ing and control scheme by which a deformed HALO optical system is corrected in 
near real time. In one of the techniques worked on by the Itek Corporation, 
the wavefront errors of the deformed system are measured with a wavefront sen¬ 
sor for several object points near the edge of the field-of-view of the sys¬ 
tem. The wavefront sensor is placed around the focal plane. The wavefront 
errors are decomposed into mirror figure errors and actuator displacements are 
generated to correct them with a deconvolution algorithm. This report de¬ 
scribes a limited blind test that CSDL prepared for Itek to test and validate 
this algorithm. We conclude that the algorithm successfuly determined the 
actuator displacements from the wavefront errors considered in this test. The 
residual wavefront errors were negligibly small and below the wavefront sensor 
noise. 


It is recommended that one more blind test be prepared for Itek's 
deconvolution algorithm to include more realism. Some of the desirable fea¬ 
tures of this test are listed below. 


1) Mirror deformations will not be produced by tho actuators used for 
correction, because, in practice, a mirror is not deformed by the 
movement of actuators. 


2) Some static aberrations will be added to the optical system to 
simulate a nondiffraction-limited undeformed system. 







The on-axis wavefront error data may not be excluded, since the 
wavefront sensor placed around the HALO focal plane cannot sense 
it. 

Tfte actuator influence functions used for mirror figure correction 
will be slightly different from those used by Itek in their algo¬ 
rithm, to simulate errors in the knowledge of influence function. 

The influence function for edge actuators may be different from 
that for interior actuators. 




SECTION 2 


SIMULATION EXTENSIONS PROJECT 


2.1 Introduction 


The overall objective of this project is to identify and quantify those 
generic mission, scenario and sensor parameters (as well as their interac¬ 
tions) which drive the performance of space-based staring IR surve?llance sys¬ 
tems. The approach taken to achieve mis overall objective has proceeded 
along two main avenues: Simulation Enhancements; and Simulation Applications. 

Efforts in the Simulation Enhancements area are directed toward general¬ 
izing those models in the Draper Integrated Simulations (Dis) which are too 
limited in scope or perhaps restricted to one particular system. Figure 2-1 
presents a simplified block diagram of the DIS. In support of the Simulation 
Enhancements effort, work has been performed (or is in progress) c-n models for 
the platform, sensor (including focal plane and signal processor), scene, and 
techniques for interfacing the DIS with an in-house image processor for 
graphics display. 

Utilizing the tools developed as a result of the Simulation Enhancements 
work, the Simulation Applications effort studies issues related to generic 
surveillance system performance drivers. Typical examples are sensor line-of- 
sight stability, focal D , and critical scene/scenario characteristics such 
as spectral interval, orbit, scene time-of-day, etc. 

2.2 Simulation Enhancements 


2.2.1 Platform Simulation 

Development of the Integrated Large Space Structures Simulation (ILS 3 ) 
was begun in 1977 in response to a DARPA need for simulating the optical 
performance of the HALO system. While intended to be general purpose, the 
internal structure of the ILS 3 (i.e., input, coordinate systems, interpola¬ 
tion methods, contact systems, etc.) was rigidly based on the HALO design. 

More recently, many analysis efforts have concentrated on single reflector or 
feedreflector systems. In the absence of a suitable raytrace program to pro¬ 
vide the necessary inputs for utilizing the ILS 3 , alternate special-purpose 
software packages were developed which did not require a raytrace and which 
could compute wavefront errors for these simple paraboloidal or spherical sys¬ 
tems. While these software packages were adequate for their intended pur¬ 
poses, this approach has three major disadvantages: the packages are problem- 
specific and must be modified, sometimes extensively, for other applications; 
no linkage with the DIS '<? provided; and documentation tends to be scattered 
and incomplete. 

The purpose of the enhancements to the ILS 3 is to assemble a package 
which can be used to analyze; both multi-mirror raytraced systems and single 
surface, paraboloidal or eoherical systems. This 't being done in five steps: 
incorporating many of th j previously proposed ILS ' changes; incorporating a 
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general purpose version of the single mirror software; adding a compatible 
raytrace program; simplifying data communication; and providing linkage to 
the NASTPAll and DISCOS programs. The use of NASTRAN is important both to pro¬ 
vide input to the top level ILS^ modules (dynamics, steady-state, and quasi¬ 
static) and in many cases to serve as a replacement for them. These enhance¬ 
ments will provide significant new capabilities for dealing wiJi both infrared 
and radar systems. 

An initial plan nas been formulated and preliminary analysis completed 
to accomplish the abo'/a objectives. The remainder of this effort is cur¬ 
rently in process. 

2.2.2 Sensor Simulation 


Many aspects of the focal plane and signal processor models originally 
in the DIS were specific to the Mini-HALO system. .In order to be able to 
study effectively the capabilities of a generic space-based staring infrared 
surveillance system, it was necessary to generalize certain aspects of each 
of these models. 

The principal enhancement to the focal plane model was to implement a 
menu of user-selectable noise models. In addition to the standard noise model 
corresponding to the basic Mini-HALO focal plane (including saturation and 
reset features), the user may select that of an ideal focal plane (i.e., 
scene shot noise only) or may specify either D* or D /Dg^p to characterize 
focal plane noise. In the latter case the user may further specify either 
a white noise or 1/f noise spectrum. 

The initial signal processor model in the DIS corresponded to the HALO 
Signal Processor (HSP). To increase flexibility and modularity, the software 
architecture for the Signal Processor has been restructured, leaving intact, 
however, the core HSP algorithm chain. In order to be able to deal with the 
needs of a general class of space-based staring IR surveillance systems, the 
character of many of these algorithms has been parameterized, and a number of 
new algorithms have been added. The focus of much of this activity has been 
in the following areas; target signature templates; thresholding procedures; 
threshold level computation and scaling; and "system track" and "acquired 
track" criteria. 

The Sensor Simulation Enhancements are complete at this point and now 
being used in DIS applications. 

2.2.3 Scene Simulation 


The single most important initiative undertaken in the Simulation 
Extensions Project has been the development of an in-house generic scene simu¬ 
lation in order to bring within the scope of the DIS the capability to gener¬ 
ate, manipulate and analyze terrestrial scene data sets as functions of the 
major surveillance system, scenario and mission parameters. This effort has 
been accomplished by working in conjunction with Photon Research Associates 
(PRA). Under subcontract to CSDL, PRA is developing a software simulation 
called GENESSIS to meet DIS needs for scene generation and manipulation. The 
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PRA work is proceeding in two phases: GENESSIS-I, to provide a near-term 
capability with certain simplifying assumptions and a limited scene data base; 
and GENESSIS-II, in which all critical scene parameters will be modelled and 
automated access provided to the DMA/Landsat data base to enable generation 
of scenes corresponding to any region for which DMA (Defense Mapping 
Agency) data exists# Tables 2-1 and 2-2 provide a summary of the principal 
features and current status of GENESSIS-I and GENESSIS-II, respectively. 

Figure 2-2 presents a simplified block diagram of the GENESSIS software 
architecture. 

2.3 Simulation Applications 

The Simulation Enhancementsdescribed above have considerably expanded 
the scope and capability of the DIS *"o conduct end-to-end performance assess¬ 
ments of space-ba? vd staring infrared surveillance systems. As an illustra¬ 
tion of this nev apability, the DIS has been used to study the effect of 
choice of spectral interval and selection of scene local time of day on clut¬ 
ter leakage through the signal processor. The basic approach was to begin 
with the California coast (Santa Cruz) data base and to use GENESSIS to pro¬ 
duce two scenes in the 3.6-4.0 pm band and two scenes in the 8.0-9.0 pm band, 
each pair of scenes corresponding to 8 AM and to noon local time, respec¬ 
tively. A nominal surveillance system mission and configuration were assumed, 
and the above scenes were processed through the DIS under varying levels of 
line-of-sight (LOS) jitter. Tables 2-3 and 2-4 summarize the scenario, sensor 
and signal processor configurations for this investigation. These parameter 
values reflect the basic configuration for a system whose principal objective 
is the detection of strategic aircraft against a terrestrial background. The 
figure of merit used to measure clutter leakage in the signal processor was 
average number of threshold exceedances per frame. Figure 2-3 illustrates 
the four scenes used. The spectrum at the bottom of the figure indicates 
relative temperature (or scene radiance) within each scene. 

Since the average number of threshold exceedances per frame is a func¬ 
tion of threshold level in the signal processor as well as the level of LOS 
jitter, the results of this study may be presented from those two points of 
view. Figure 2-4 shows the dependence on LOS jitter level of the average num¬ 
ber of threshold exceedances per frame. The particular threshold level 
selected is representative of what would be used to acquire a target with a 
relatively high signal-to-noise ratio. The general trend of the data is such 
that, for a given level of LOS jitter, higher numbers of exceedances are found 
in the LWIR band (i.e., 8.0-9.0 pm) than in the SWIR band fi.e., 3.6-4.0 pm). 
In either band, the higher number of exceedances always occurs at noon local 
time as opposed to 8 AM. 

Figures 2-5 through 2-7 show the dependence of average number of thresh¬ 
old exceedances per frame on threshold level, for progressively increasing 
values of LOS jitter. Many of the same trends are also evident in this data. 
Decreasing the threshold results in a rapid increase in threshold exceedances. 
At a given threshold value, higher numbers of threshold exceedances always 
occur in the LWIR band as opposed to the SWIR, and, within each band higher 
numbers of exceedances occur at noon in comparison with 8 AM. Increasing the 
LOS jitter level results in a translation upward of the curves at higher 
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Table 2-1. Generic scene simulation Phase I (GENESSIS-I) 


FEATURES 

• User-specified variables 

-Observer altitude, zenith angle 
-Spectral interval (2.5-13.0 ym) 

-Atmospheric model (LOWTRAN) 

-Field-of-view location (within overall scene) 

-Scene angular resolution 

• Five representative terrain data bases 

-California Coast (near Santa Cruz) 

-Brooks Range Mountains of Alaska 
-Arctic Tundra 
-Middle East 
-Central Europe 

• Each data base topology extracted from DMA data base 

• Each data base registered with Landsat data to obtain 
materials assignments 

-Limited to 14 material types 

• Two representative cloud templates (for superposition over 
terrain data bases) 

• Scene size: 40 km x 40 km 


STATUS 


Operational at CSDL 

















Table 2-2. Generic scene simulation Phase II (GENESSIS-II) 


FEATURES 

• Enhanced data base 

-Material types expanded to 20 

-Geometrical representations for three additional cloud 
patterns (Derived from NOAA data) 

• GENESSIS/DMft'-Landsat interface 

-Software package to automate generation of input data bases 
from raw DMA and Landsat data 

-Will permit generation of any scenes for which DMA data 
exists 

• Generalized treatment of key atmospheric effects 

-Continuously variable surface-level parameters (e.g., 
temperature, humidity, wind velocity, skyshine, solar 
scattering) 

• Optimized software architecture 

• Detailed user's manual 


STATUS 

• GENESSIS-II (Part 1) subcontract to PRA currently underway 
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Table 2-3. Scenario and sensor specification for GENESSIS 
Santa Cruz investigation. 


PARAMETER 

VALUE 

Target 

None 

Backgrounds 


-Location 

Santa Cruz (GENESSIS) 

-Spectral Bands 

3.6-4.0 pm 

8.0-9.0 pm 

-Times of Day (Local) 

8 AM 

NOON 

Spacecraft 


-Drift 

None 

-Jitter 

0.1 prad rms 

0.5 prad rms 

1 .0 prad rms 

-Altitude 

35,700 km 

-Boresight 

NADIR 

Optics 


-Aperture Diameter 

1.4 m 

-Obscuration Ratio 

0.35 

-Thruput 

0.135 

-Point-spread Function 

Gaussian (c = 5.01 prad) 

Focal Plane 


-Detector Geometry 

25 x 25 Array of 100 pm 


(20 prad) Square Pixels 

-Frame Duration 

3.4 sec 

—D* 

2.8 (10) 12 cm /Hz/W 

-Responsivity non-uniformity 


.Fixed Pattern 

0% 

.Random 

15% 


0 




















Table 2-4. Signal processor configuration for GENESSIS 
Santa Cruz investigation. 


PARAMETER 

VALUE 

Difference Filter 

3rd Order 

Threshold Level 

User Input 

Assumed Target Contrast 

Negative 

Threshold Templates 

CSDL MTD 

Threshold Algorithm 

CSDL Rationalized 
Exceedance Thresh¬ 
olding (Version 1) 

Cluster Basis 

Threshold Exceedances 

Centroid Basis 

Threshold Exceedances 

System Track Criterion 

3 out of 3 

Acquired Track Criterion 

10 out of 12 

System Track 

Threshold Scaling 

0.5 

Least Squares Track 
Linearity Criterion 

±1 .8 Pixels 





































AVERAGE NUMBER OF EXCEEDANCES PER FRAME 



Figure 2-4. Average clutter leakage as a function of line-of-sight 
jitter level for the GENE3SIS Santa Cruz scene, with 
local time-of-day and spectral band as parameters. 







Figure 2-5, Average clutter leakage as a function of threshold level 
for the GSrJESSIS Santa Cruz scene, with local time of day 
and spectral band as parameters. (LOS Jitter =0.1 prad 















iH 



| X - 8.0 - 9.0 pm 


LOS JITTER 
FRAME SIZE 

NUMBER OF DATA FRAMES 
FRAME DURATION 


= 0.5 prud rms 
= 25 X 25 PIXELS 
= 25 
= 3.4 s 


X = 3.6 - 4.0 pm 


THRESHOLD VALUE 


Figure 2-6. Average clutter leakage as a function of threshold level 
for the GENESSIS Santa Cruz scene with local time of day 
and spectral band as parameters. (LOS Jitter =0.5 prad 
RMS) 







AVERAGE NUMBER OF 
EXCEEDANCES PER FRAME 


X ■ 8.0 - 9.0 



THRESHOLD VALUE 


Figure 2-7. Average clutter leakage as a function of threshold level 
for the GENESSIS Santa Cruz scene, with local time of day 
and spectral band as parameters. (LOS Jitter = 1.0 lirad 
RMS) 


16 





threshold values, as can be seen by comparing Figures 2-5 and 2-6* In addi¬ 
tion, there is a compression within each spectral band of the exceedance 
curves corresponding to noon and to 8 AM* In Figure 2-7, an LOS jitter level 
of 1.0 yrad rms (ten times the value used in Figure 2-5) results in the co¬ 
alescing of these two curves for the LWIR band. 

A subject of great interest in the study of signal processor performance 
and scene/sensor interaction is a determination of the conditions under which 
false target tracks arise. Figure 2-8 summarizes graphically the false tracks 
which resulted during the course of the present study. Separate numbers, 
i.e., "1" versus "2”, are used to designate the relative position of distinct 
tracks within the overall 25 x 25 pixel field-of-view. Consistent with the 
data presented in Figures 2-4 thi »ugh 2-7, false tracks were found at low 
threshold levels where the average number of exceedances is very high. Under 
these conditions the likelihood of a time sequence of individual exceedances 
occurring in such a fashion as to satisfy the signal processor track formation 
and acquisition criteria is considerably increased. 

2.4 Conclusions 

The focal plane and signal processor simulation enhancements are com¬ 
plete, and the platform simulation enhancements are in process and proceeding 
satisfactorily. Phase I of the Generic Scene Simulation (GENESSIS) activity 
is complete and operational at CSDL. The last phase of this activity is cur¬ 
rently in process under subcontract to CSDL at Photon Research Associates. 

The Draper Integrated Simulations thus represent a highly sophisticated tool 
for end-to-end performance trade-offs of a variety of space-based systems, and 
especially staring mosaic IR surveillance systems. 

LOS-induced clutter leakage through the signal processor has been simu¬ 
lated as a function of spectral interval and scene time-of-day for a generic 
staring mosaic IR surveillance system. The functional dependence of the 
scene/sensor interaction is complex and has been described quantitatively by a 
series of curves showing average number of threshold exceedances per frame 
both as a function of LOS jitter for a particular threshold and as a function 
of threshold for particular LOS jitter values. Qualitatively it is clear that 
increased LOS jitter yields increased clutter leakage in all cases consid¬ 
ered. Clutter leakage in the short-wave infrared (SWIR) band of 3.6-4.0 ym is 

significantly less than in the long-wave infrared (LWIR) band of 8.0-9.0 ym. 

For either SWIR or LWIR, clutter leakage at noon (local time) is greater than 

early morning. For a fixed threshold value, increased LOS jitter drives early 
morning clutter levels toward noon clutter levels in both the SWIR and the 
LWIR bands. And, finally, low threshold values in the presence of LOS jitter 
can result in the generation of false acquired tracks. 
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SECTION 3 



HALO OPTICS 
DECONVOLUTION TEST 


3.1 Introduction 


HALO is a multi-mirror large optical imaging system with a wide field 
of view. To maintain alignment as well ac. image quality, the mirrors are 
actively controlled with actuators. The teim "dec...:volution" implies a 
sensing and control scheme by which a deformed HALO optical system i-~ cor¬ 
rected in near real time. In one of the techn ques worked on by the Itek 
Corporation, the wavefront errors of the deformed system are measured with a 
Wcvefront sensor for several object points near the edge of the field-of- 
view of the system. The wavefront sensor is placed around the focal plaae. 
The wavefront errors are decomposed into mirror figure errors and actuator 
displacements are generated to correct them with a deconvolution algorithm. 
This report describes a limited blind test that CSDL prepared for itek to 
test and validate this algorithm. We conclude that the algorithm success¬ 
fully determined the actuator displacements from the wavefront errors. 

3.2 Optical System for Deconvolution Test 

For a hardware demonstration of their deconvolution algorithm, Itek 
has been planning to use an optical imaging system consisting of three 
active flat mirrors and a pair of lens doublets mounted back-to-back in a 
1:1 magnification. A schematic of this system is shown in Figure 3-1. 

The lenses are the imaging elements of the system, and the mirrors 
are the deformable elements. Aberrations are introduced into the system by 
deforming the mirrors. The active mirrors, thus simulate the deformable 
aspects of a HALO system. The footprints on the mirrors of rays from a 
given point object are similar to those on the S-HALO or WALRUS mirrors. A 
perspective view o.f the optical system is given in Figure 3-2. It shows the 
coordinate system chosen to designate the actuator and ray distributions. 

3.3 Actuator Distribution and Influence Function 


Each of the three mirrors has 37 actuators distributed on an equi¬ 
lateral triangular lattice as shown in Figure 3-3. In consultation with 
Itek, the influence function of the actuators was chosen to be similar to 
the one for active mirrors developed for the compensated imaging program. 
This influence function is given by 


f(r) = [l - 1.2(r/d) 2 ] 


r < d//i72 


» 0 , r > d/ZTTT , 


/ 


where 
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r = distance from the actuator, 
d = actuator spacing. 

The variation of f(r) with r is shown in Figure 3>-4. It is noted 
that the influence function is radially symmetrical about the actuator 
position. Moreover, one actuator does not produce any influence on 
another. The factor of 1.2 was chosen to minimize the effect of quilting 
when all actuators are displaced by the same amount. With a 1.2 factor, 
quilting is at most 8 %. When an actuator is displaced by an amount h, the 
deformation produced in the mirror is given by hf(r). In our problem d = 
3cm. 


3.4 Mirror Deformations 


The deformations of the mirrors used in the deconvolution test are 
given in Figure 3-5. The numbers in this figure give the actuator 
displacements in units of X/100. Mirror M^ was deformed according to 
cryogenic deformations of an Itek HALO mirror. Mirror M 2 has deformations 
described by Zernike polynomials representing coma and astigmatism. Mirror 
M 3 was given a rigid-body displacement along its normal. The mean, 
peak-to-peak, and standard deviation of the actuator displacements are 
indicated in the figure for each mirror. 

3.5 Selection of Point Objects 

Nine point objects were selected for ray tracing and determining the 
system wavefront errors. Eight of these are uniformly distributed around 
the edge of the field-of-view of the system (1.5° from the axis), and one is 
placed on the system axis as indicated in Figure 3-6. The on-axis point 
object was used only for the sake of completeness to include the deforma¬ 
tions of the central portion of mirror Mj. The wavefront sensors placed 
around the focal plane cannot sense the effect of these deformations. In 
practice, the central portion of M^ will, for example, be corrected by 
CSDL's image sharpening approach. ' 

3.6 Mirror Footprints 

Figure 3-7 shows *'he footprints on the mirrors of the image forming 
rays from point objects 1 and 2. We note that, as in a HALO optical system, 
the footprints of rays from point objects such as 2 , near the edge of the 
field of view, on mirrors Mj and M£ cover the mirror portions used in imag¬ 
ing. However, their footprint on mirror Mj is not only small compared to 
the mirror size, but also it is off the mirror center. Accordingly, 
the footprints of the rays from the 8 off-axis point objects cover the outer 
portion of the mirror M} but not its central portion. 

The image of the on-axis point object is formed at the center of the 
focal plane. The central portion of mirror M^ can be corrected by the image 
sharpening technique using the focal plane data. 

































Figure 3-5b. Mirror 
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3.7 


Wavefront Errors 


When the optical system deforms, the wavefront error introduced at a 
certain point on the wavefront is given by the change in' the path length of 
a ray passing through that point. The change in the path length of a ray 
from a given point object in terms of the mirror deformations is given 
according to 3 


AWj. = 


3 

I 

m=1 


n. 

lm 


Sr. 

im 


cos 0 


im' 


where 

A 

n. = unit vector along the normal to the surface of undeformed mirror 
m at the point of incidence of the ray i, 

5r^ m 3 displacement of the point of incidence of ray i on mirror m, 

0 ^„ n 3 angle of incide ice of ray i on the undeformed mirror m. 

Since we are dealing with flat mirrors, ni m for a given mirror is indepen¬ 
dent of i. 

Figure 3-8 shows the distribution of 89 rays in the pupil used for 
determining the wavefront errors. The wavefront errors corresponding to 
point objects 1 and 2 when the optical system is deformed according to 
mirror deformations of Figure 3-5 are given in Figure 3-9. These wavefront 
errors include X/20 random error to simulate the noise of a wavefront 
sensor. We provided Itek .with noisy wavefront errors at 89 points for each 
of the nine point objects. 

3.8 Test Results 


3.8.1 Mirror Deformations 


From the wavefront error data for the nine point objects we provided 
to Itek, they determined the actuator corrections (displacements) required 
:o correct the figure of the mirrors to reduce the wavefront errors. The 
actuator corrections determined by their deconvo.lution algorithm are given . 
in Figure 3-10. The residual errors obtained by subtracting these correc¬ 
tions from the errors given in Figure 3-5 are also shown in Figure 3-10. 

The mean and the standard deviation of the residual displacements are listed 
in the figure for each of the three mirrors. We note that the residual 
errors are quite small except for the corner actuators. The region of 
mirrors where corner actuators are located is not used for imaging. Hence, 
the effect of these actuators on the mirror figure is not observed by the 
wavefront sensor. 
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3.8.2 Wavefront Errors 


The wavefront errors of the corrected system were determined by 
tracing 89 rays as in the case of the deformed system. These errors were 
quite small (c^ < X/40) indicating an excellent correction. Random noise 
of X/20 was added to these errors to simulate the wavefront errors that a 
wavefront sensor.would measure. The noisy residual wavefront errors for 
point objects 1 and 2 are shown in Figure 3-11. The mean, peak-t.o-peak, and 
standard deviation of the wavefront errors are also given. It is evident 
from the standard deviation that most of the error is indeed due to the 
wavefront sensor. 

3.8.3 Point-Spread Functions 


The aberrated and the corrected point-spread functions for point 
objects 1 and 2 are shown in Figures 3-12 and 3-13, respectively. The 
standard deviation of the aberration, the Strehl ratio (aberrated irradiance 
at the center relative to its unaberrated value) and peak irradiance are 
also given in the figures. For example, for point object 1, the deformed 
system gives a Strehl ratio of only 0.025. The light is spread over a region 
31XF wide compared with approximately 2XF for the corrected (or aberration- 
free) system. With a correction determined by Itek deconvolution algorithm, 
a Strehl ratio of approximately 0.95 is obtained. The corrected Strehl 
ratios in Figures 3-12 and 3-13 are slightly lower because of the wavefront 
sensor noise. 

3.9 Discussion, Conclusions, and Recommendations 


An optical system using three flat mirrors as deformable elements has 
been used for testing the deconvolution algorithm developed by Itek. Each 
mirror has 37 actuators, and 89 rays have been used to determine wavefront 
errors of a deformed system. CSDL provided Itek with a known-deformation 
test and a blind test. The known-deformation test was prepared so that CSDL 
and Itek softwares could be tested for compatibility. On the blind test, 
two of the 9 wavefronts that Itek was given are shown in Figure 3-9. 

Itsk's deconvolution algorithm determined the actuator displacements 
on each mirror and these were subtracted from the input displacements given 
in Figure 3-5. The corrected system was ray-traced to determine its imaging 
quality. It gave a Strehl ratio of approximately 0.9 indicating an 
excellent correction. 

We conclude that Itek's deconvolution algorithm successfully 
determined the actuator displacements from the wavefront errors. The 
residual wavefront errors were negligibly small and below the wavefront 
sensor noise. 

We recommend that one more blind test be prepared for Itek's decon¬ 
volution algorithm to include more realism. Some of the desirable features 
of this test are listed below. 














































are the optical 



















































1. Mirror deformations will not be produced by the actuators used 
for correction, because, in practice, a mirror is not deformed by 
the movement of actuators. 

• . % 

2. Some static aberrations will be added to the optical system to 
simulate a nondiffraction-limited undeformed system. 

3. The on-axis wavefront error data may not be provided, since the 
wavefront sensor placed around the HALO focal plane can not sense 
it. 

4. The actuator influence functions used for mirror figure 
correction will be slightly different from those used by Itek in 
their algorithm, to simulate errors in the knowledge of influence 
function. 

5. The influence function for edge actuators may be different from 
that for interior actuators. 

6 . A segmented mirror will be used. 

7. The mirror size will be approximately the same as the illuminated 
region. 
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